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摘要 :基于 1990 一 2013 年 土地 利用 数据 和 统计 年 鉴 数据 ,构建 CPI 指数 修订 生态 系统 服务 价值 (ESV ) 系数 ,并 进行 准 河上 游 士 
地 利用 结构 及 空间 格局 的 时 空 演变 特征 分 析 和 ESV 评价 研究 ,探讨 ESV 对 土地 利用 结构 和 格局 变化 的 响应 关系 。 结 果 表 明 : 
(1) 淮 河上 游 区 域 ESV 整体 呈 上 升 趋势 ,1995 一 2005 年 ESV 出 现 波动 下 降 , 在 2010 年 前 后 明显 上 升 ,2010 一 2013 年 ESV 增值 
高 于 其 他 时 段 ;(2)ESV 增 量 呈 现 山地 > 丘陵 区 > 平原 区 的 地 形 梯度 效应 , 且 其 生态 增值 主要 来 源 于 林地 和 水 域 贡 献 ;(3) ESYV 与 
土地 利用 时 空 变化 和 土地 景观 格局 丰 度 和 聚集 度 变 化 存在 明显 相关 ,说 明了 合理 的 土地 利用 开发 方式 及 适当 的 开发 速度 将 有 
助 于 改善 生态 系统 服务 ,开发 过 程 要 兼顾 地 类 优势 ,避免 地 类 斑 块 过 度 破碎 化 ;(4) CPI 指数 修订 ESV 系数 方法 ,能 够 消除 时 间 
序列 内 物价 变动 因素 ,使 得 计算 结果 更 贴近 实际 ESV ,提高 了 多 个 时 期 生态 系统 服务 评价 的 准确 性 和 可 比 性 。 
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Abstract: It is important to study the effect of land use change on ecosystem service value (ESV). Most previous studies 
have tended to directly use the ecosystem service equivalent value to calculate the ESV and then have compared multi period 
ESVs. Because of the use of different price bases, even if the value of the same ecosystem is assessed using multi period 
ESVs, the conclusions can be very different, which means there has to be a deviation evaluation. In order to improve the 
accuracy of these comparisons, we developed a Consumer price index( CPI) revised ESV coefficient method for multi period 
ESV computation. The upper reaches of the Huaihe River was used as a case study to estimate the spatial distribution and 
dynamics of ecosystem services. In this study, we used the land use and statistical yearbook data to investigate how the 
spatial-temporal evolution of land use structure and spatial patterns affect the ESV. The results showed that (1) the ESV 
generally showed an increasing trend in the study region. Between 1995 and 2005, the ESV decreased, but then increased 
significantly around 2010. Furthermore, the increase from 2010 to 2013 was higher than that for any other period. (2) The 
ESV increase was higher in mountain areas than in hilly areas, and the increase in hilly areas was higher than that in the 


plains, which shows a topographical gradient effect. The increase was mainly because of the contributions made by forest and 
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water ecosystems. (3) There were significant correlations between the ESV and land use change and Aggregation index 
(AI) and Number of patches( NP ) landscape pattern indexes. These correlations show that a reasonable model for land use 
development and an appropriate development rate will help enhance ecosystem services. The development process should 
take into account the advantages that different land types can offer so that excessive fragmentation of land type patches can 
be avoided. (4) A method that uses the CPI index to revise ESV coefficients means that the values more closely reflect the 


actual ESV and improve multi-period ecosystem service evaluation comparability and accuracy. 


Key Words: ecosystem services value; land use change; landscape pattern; spatial-temporal heterogeneity; the upper 


reaches of the Huaihe River 


在 经 济 高 速 发 展 的 影响 下 ,土地 资源 的 交易 和 占用 愈 发 活路 ,促使 土地 使 用 类 型 和 利用 方式 快速 转 

变 ''"9 ,对 生态 系统 服务 价值 ( ESV， Ecosystem Services Value) 产生 了 显著 影 响 ' 引 。 自 1997 年 Costanza 在 

《Nature》 提 出 ESV 估算 原理 及 方法 后 ”… ,我 国学 者 开始 寻找 更 适合 于 中 国 ESV 评估 的 方法 。 比 较 有 代表 性 
本 研究 有 :欧阳 志 云 等 引 综 合 运用 生态 学 及 经 济 学 方法 ,探讨 了 大 区 域 生态 服务 功能 的 内 涵 与 评价 方法 。 陈 仲 
SU ”新 等 (1 \ 毕 晓 丽 等 "1 采用 Costanza 等 人 提出 的 生态 系统 类 型 单位 价值 对 中 国生 态 系统 功能 与 效益 进行 了 价值 
估算 。 何 浩 等 !5 结合 遥感 技术 计算 了 2000 年 中 国 陆地 ESV。 谢 高 地 等 '| 参考 Costanza 等 人 部 分 成 果 , 同 时 
综合 中 国 专业 人 士 的 生态 问卷 调查 结果 ,建立 了 中 国 陆地 生态 系统 单位 面积 服务 价值 表 , 在 此 基础 上 ,于 
2006 年 对 中 国 500 位 具有 生态 学 背景 的 专业 人 员 进 行 问卷 调查 ,得 出 新 的 生态 系统 服务 评估 单价 体系 ,并 
于 2015 年 进行 修订 "1 , 构 建 了 年 内 逐 月 和 空间 上 分 省 的 时 空 动 态 当量 因子 表 。 虽 然后 续 也 有 不 少 学 者 进行 
ESV 方法 的 改进 与 应 用 研究 于 2 ,其 实质 仍 以 谢 高 地 等 于 2003 年 提出 的 基于 单位 面积 价值 当量 表 为 基础 
进行 。 
GAN 由 于 计算 ESV 的 价格 基础 不 同 ,出 现 基 于 单位 面积 ESV 当量 的 计算 结果 远 远 高 于 基于 单位 面积 价值 量 
~ 的 计算 结果 "中 ,直接 弱化 ESV 在 土地 管理 实践 中 的 应 用 '“] 。 张 舟 等 "5 通过 分 析 Zang 等 "两 种 ESV 当 
量 下 的 ESV 变化 情况 ,发 现 ESV 对 当量 参数 敏感 性 小 的 特点 ,并 认为 边际 变化 方法 可 用 于 ESV 变化 研究 。 
因此 , 现 有 文献 常 对 两 个 时 期 ESV 做 差 值 ,进而 讨论 ESV 的 时 序 演变 情况 !*221。 事 实 上 ,在 排除 各 ESV 当 
.三 量 系 数 影响 敏感 性 小 情况 下 ,单位 面积 价值 成 为 引起 ESV 变化 的 主要 因子 。 现 有 单位 面积 价值 用 研究 区 主 
要 粮食 单位 面积 产值 的 1/7 进行 蔡 代 ,该 部 分 计算 并 未 考虑 物价 上 涨 因素 ,造成 后 续 进 行 边际 变化 分 析 的 
ESV 时 序 差 值 仍 不 是 同一 定 基 ,从 而 导致 即使 生态 系统 实际 价值 相同 也 会 表现 出 以 评估 年 份 价格 为 表征 的 价 
值 量 不 同 现象 "| ,直接 影响 边际 变化 分 析 结果 准确 性 。 另 外 ,已 有 研究 基于 土地 面积 进行 ESV 变化 研究 ， 
更 多 关注 于 土地 资源 面积 变化 问题 ,而 忽略 土地 利用 面积 变化 所 引发 的 结构 与 格局 的 变化 中 ,特别 是 土地 
利用 过 程 变 化 。 即 使 一 些 学 者 研究 农 地 非 农 化 的 代价 性 损失 站, 仍 未 对 土地 利用 时 空 变化 过 程 与 ESV 变化 
的 影响 关系 进行 分 析 与 探讨 。 
淮河 流域 地 处 我 国 南北 气候 过 渡 带 和 中 原 腹 地 。 至 2010 年 ,流域 以 不 足 全 国 13% 的 耕地 面积 和 19% 的 
播种 面积 ,承载 了 全 国 20.4% 的 乡村 人 口 ,生产 了 全 国 23% 的 粮食 ,提供 了 全 国 14% 的 农业 产值 '”。 在 此 自 
然 资源 敏感 性 和 经 济 需求 迫切 性 并 存 的 背景 下 ,本 研究 选取 淮河 上 游 区 域 ,进行 土地 利用 及 景观 格局 演变 对 
ESV 影响 机 理 研究 ,提出 基于 单位 面积 生态 系统 服务 价值 CPI( 居民 消费 价格 指数 ) 系数 修订 方法 ,并 结合 
ESV 损益 方法 ,进行 ESV 过 程 变 化 研究 。 则 在 从 宏观 和 微观 空间 尺度 上 ,定量 评估 土地 利用 变化 引发 的 生态 
效应 ,以 期 为 淮河 上 游 未 来 土地 资源 可 持续 利用 和 社会 经 济 可 持续 发 展 提供 理论 和 方法 支持 。 


1 研究 区 概况 及 研究 方法 


1.1 研究 区 概况 
淮河 上 游 处 于 我 国 南 北 气候 的 过 渡 带 。 在 气候 区 划 中 ,以 淮河 和 苏 北 灌溉 总 渠 为 界 , 北 部 属 暧 温带 半 湿 
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润 地 区 ,南部 属 亚热带 湿润 区 ,年 平均 气温 为 11 一 
16%C ,并 呈现 出 自 南 向 北 逐 渐 递 减 趋势 。 流 域内 降水 时 
空 分 布 特征 明显 ,具有 南 多 北 少 ,山区 高 于 平原 ,年 内 降 
水 多 发 生 于 汛期 ,年 际 降水 量变 化 大 的 特点 。 淮 河上 游 
南部 由 西向 东 依 次 为 桐 柏 山 和 大 别 山 ,山地 丘陵 占 优 
势 , 海 拨 1000m 以 上 的 中 山 面积 约 2000km ,海拔 
1000m 以 下 的 低 山 丘陵 面积 近 30000kmz 2 。 淮 河上 游 
北部 主要 为 平原 ,是 我 国 重要 的 粮食 .棉花 .油料 生产 基 S| 
地 ,农业 生产 水 平和 土地 利用 率 较 高 。 2 

淮河 上 游 流 域 包括 信阳 市 .平顶山 市 .南阳 市 和 驻 
马 店 市 部 分 地 区 ,为 便于 社会 经 济 资料 的 获取 分 析 及 
成 果 的 应 用 ,本 文 把 政 区 面积 60% 以 上 在 淮河 上 游 流 ”31。 
域 边界 内 的 县 (市 .区 ) 纳 入 研究 范围 。 
1.2 ”数据 来 源 

研究 所 用 流域 行政 区 划 数 据 和 6 个 年 份 的 土地 利 
用 数据 , 均 来 自 于 国家 科技 基础 条 件 平 台 一 一 国家 地 球 
系统 科学 数据 共享 平台 一 黄河 下 游 科学 数据 中 心 ,30m DEM 数据 来 源 于 中 国 科 学 院 数 据 云 平台 。 计 算 各 
生态 系统 单位 面积 服务 价值 所 用 耕地 数据 粮食 产值 数据 来 源 于 1990 年 .1995 年 .2000 年 .2005 年 .2010 年 
和 2013 年 河南 统计 年 鉴 。 消 除 物价 上 涨 所 用 的 CPI 指数 数据 ,来源 于 国家 统计 局 网 站 。 土 地 利用 数据 源 为 
landsatTM 影像 ,采用 人 机 交互 解 译 方法 获得 ,数据 准确 率 达 到 90% ,满足 研究 需要 。 
1.3 数据 处 理 方 法 

土地 利用 数据 主要 借助 于 ENVI 5.1 和 ArcGIS 10.2 平台 进行 影像 处 理 和 空间 分 析 计 算 ,景观 指数 采用 
Fragstats 4.2 软件 进行 获取 ,参与 计算 的 栅 格 数据 单元 为 30mx30m。 相 关系 数 计算 基于 SPSS 统计 软件 。 土 
地 利用 变化 分 析 所 用 到 的 变化 率 转移 矩阵 计算 公式 参照 文献 17 ,景观 格局 指数 计算 方法 详细 说 明 参 照 文献 
22。 下 面 着 重 介绍 生态 服务 价值 修订 及 计算 方法 。 

(1) 生 态 服 务 价值 


图 1 研究 区 位 置 示意 图 


Fig.l Sketch map in the study area 


ESV = > A,VC, (1) 
k=1 


式 中 ,ESV 为 研究 区 生态 系统 服务 总 价值 (元 ) ;n 为 土地 利用 类 型 的 数量 ;4 为 研究 区 内 土地 利用 类 型 有 的 分 
布 面积 (km ) ;VC 为 第 类 土地 利用 类 型 的 ESV 系数 , 即 单位 面积 上 土地 利用 类 型 大 的 ESV( 元 km ”a )。 
在 式 (1) 的 基础 上 ,得 出 ESV 损益 计算 公式 : 
Py = (VC; ~ VC) A; (2) 
式 中 ,P; 为 第 i 类 土地 利用 类 型 转化 为 第 j 类 土地 利用 类 型 后 的 ESV 损益 ;VC;, VC, 分别 为 第 i 类 土地 利用 类 
型 的 ESV 系数 ;4 为 第 i 类 土地 利用 类 型 转化 为 第 j 类 土地 利用 类 型 面积 。 
(2) 单 位 面积 ESV 修订 


w = CPI/100 (3) 
MT = (4) 
do ai Q2 CQ 1 Qo 
OQ, 
VC, = VC,’ x 一 (5) 
Wo 
VCx' = 当量 x 单位 面积 价值 (6) 
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式 中 ,CPI 为 居民 消费 指数 ;a, 为 去 除 单位 CPI, 一 为 第 一 期 环 比 CPT, 一 为 现 期 同 比 定 基 CPI; 当量 数据 参考 
文献 表 10 ;单位 面积 价值 为 研究 区 当年 主要 粮食 作物 单位 面积 产值 的 1/7。 
2 土地 利用 结构 和 空间 格局 时 空 变化 分 析 


2.1 土地 利用 结构 的 时 间 变 化 特征 分 析 

土地 利用 变化 主要 从 土地 利用 类 型 面积 变化 和 变化 速率 两 个 方面 进行 。 由 表 1 可 以 看 出 ,23 年 间 研 究 
区 建筑 用 地 面积 比例 持续 上 升 ,耕地 和 草地 面积 呈 下 降 趋 势 。 至 2013 年 底 , 耕 地 占 研 究 区 67.2% 的 面积 , 仍 
据 主导 地 位 ,草地 所 占 比例 由 1990 年 的 3.02% 下 降 至 2013 年 2.16%。 以 每 5a 为 间隔 计算 各 时 间 段 的 土地 利 
用 综合 变化 动态 度 , 分 别 得 到 0.48% (1990 一 1995 年 )、0.36% (1995 一 2000 年 )、0.09%(2000 一 2005 年 )、 
0.38%(2005 一 2010 年 ) 和 0.11%(2010 一 2013 年 ) ,从 中 发 现 ,1995 年 以 前 是 一 个 地 类 快速 转化 的 时 期 ,1995 
年 之 后 地 类 变化 速度 减 慢 ,2005 年 之 后 逐渐 加 速 转化 。 


表 1 淮河 上 游 1990 一 2013 年 土地 利用 类 型 构成 及 变化 率 


Table 1 Land use structure and the changing rate of the upper reaches of the Huaihe River from 1990 to 2013 


土地 类 型 土地 利用 类 型 构成 Land use type structure/% 土地 利用 变化 率 Land use change rate/% 

Land use type 1990 1995 2000 2005 2010 2013 ”1990 一 1995 1995 一 2000 2000 一 2005 2005 一 2010 ”2010 一 2013 
浊 地 Framland 69.51 70.23 69.89 69.56 67.71 67.2 0.21 -0.10 -0.09 -0.53 -0.25 
林地 Forestland 16.38 17.74 16.29 16.23 17.63 17.62 1.67 -1.64 -0.07 1.72 -0.02 
草地 Grassland 3.02 1.19 2.23 2.61 2.16 2.16 -12.11 17.44 -0.64 -0.01 -0.06 
水 域 Waters 2.86 2.28 3.01 3.31 3.26 3.26 -4.07 6.36 2.01 -0.27 0 
Mr 8.23 8.55 8.59 8.74 9.24 9.77 0.79 0.09 0.37 1.13 1.91 
ee 0 0 0 0 0 0 0 0 0 0 0 


表 2 地 类 转移 矩阵 说 明了 1990 一 2013 年 各 地 类 的 互 转 情况 ,其 中 ,B% 为 初期 转 出 比重 ,C% 为 末期 转 入 
比重 。23 年 间 耕 地 转 出 面积 最 大 ,为 1801.11km , 按 转 出 面积 多 少 依 次 为 :草地 379.82km ,林地 340.61km”， 
水 域 63.51km ,建筑 用 地 31.46km ; 转 出 比重 最 大 的 为 草地 , 转 出 率 为 34.52% ,其 次 为 耕地 7.12%。 转 人 面积 
最 大 的 为 建筑 用 地 , 达 1034.93km , 转 和 人 率 为 16.66% ,其 他 地 类 依次 为 :林地 791.10km ,耕地 496.34km ,水域 
229.73km 和 草地 64.35km ”。 因 为 未 利用 地 本 身 基数 小 ,微小 的 数量 变化 都 会 导致 较 大 的 比重 差异 ,23 年 间 
只 有 耕地 有 少量 转 为 未 利用 地 外 ,未 利用 地 并 未 发 生 其 他 变化 ,因此 ,忽略 未 利用 地 的 转移 分 析 。 


表 2 淮河 上 游 1990 一 2013 年 土地 利用 转移 矩阵 
Table 2 Land use Transfer matrix of the upper reaches of the Huaihe River from 1990 to 2013 


本 地 | 转 出 合计 
土地 利用 类 型 林地 草地 水 域 建筑 用 地 未 利用 地 Wp 

Cultivated Transfer-out 1990 年 
Land use type Forest land Grass land Waters Build-up land Unused land : 

land summattion 

耕地 /km? 23513.08 641.03 17.01 180.62 962.38 0.06 1801.11 25314.19 
B% 92.89 2:53 0.07 0.71 3.80 0.01 7.12 100.00 
C% 96.08 2.62 0.07 0.74 3.93 0.01 
林地 /km? 235.07 5624.54 42.09 24.09 39.36 0.00 340.61 5965.15 
B% 3.94 94.29 0.71 0.40 0.66 0.00 $71 100.00 
C% 3.66 87.67 0.66 0.38 0.61 0.00 
草地 /km? 196.84 133.92 720.61 20.83 28.24 0.00 379.82 1100.43 
B% 17.89 12.17 65.48 1.89 2:357 0.00 34.52 100.00 
C% 25.08 17.06 91.80 2.65 3.60 0.00 
水 域 /km? 43.74 10.34 4.49 979.43 4.95 0.00 63.51 1042.94 
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耕地 i . A s 转 出 合计 
土地 利用 类 型 林地 草地 水 域 建筑 用 地 。 未 利用 地 人 王 
Cultivated Transfer-out 1990 年 
Land use type Forest land Grass land Waters Build-up land Unused land . 
land summattion 
B% 4.19 0.99 0.43 93.91 0.47 0.00 6.09 100.00 
C% 3.68 0.87 0.38 82.41 0.42 0.00 
建筑 用 地 /km? 20.69 5.83 0.75 4.18 2964.12 0.00 31.46 2995.58 
B% 0.69 0.19 0.03 0.14 98.94 0.00 1.05 100.00 
C% 0.08 0.16 0.02 0.12 83.34 0.00 
未 利用 地 /km? 0.00 0.00 0.00 0.00 0.00 0.10 0.00 0.10 
B% 0.00 0.00 0.00 0.00 0.00 100.00 0.00 100.00 
C% 0.00 0.00 0.00 0.00 0.00 62.50 
转 人 合计 /km? 
Transfer-in 496.34 791.10 64.35 229.73 1034.93 0.06 
summattion 
C% 3.92 12.33 8.20 17.59 16.66 37.50 
2013 年 /km? 24009.42 6415.64 784.96 1209.15 3999.05 0.16 36418.38 


下 ;C% :末期 转 出 比重 


B% :初期 转 出 比 习 


2.2 土地 利用 变化 的 空间 变化 特 生 


FE 分析 


上 背 助 于 ENVI 5.1 对 两 个 时 期 的 土地 利用 类 型 图 进行 土地 利用 变化 监测 制图 , 找 出 变化 区 域 , 并 在 ArcGIS 


中 按 5 种 坡度 类 别 统计 ,和 


可 LDLD 


表 3。 表 3 表明 随 坡 度 增加 ,各 地 类 转移 量 均 有 明显 下 降 , 并 呈现 地 类 自然 属性 


效应 。 其 中 ,耕地 转 林地 ,草地 转 林 地 两 类 转移 更 多 发 生 在 6 一 15° 坡 度 区 间 。 而 林地 转 耕 地 耕地 转 建筑 用 
地 则 更 多 发 生 在 0 一 2* 区 域 。 


表 3 


淮河 上 游 1990 一 2013 年 土地 转移 区 域 坡度 统计 


Table 3 Land transfer regional slope statistics of the upper reaches of the Huaihe River from 1990 to 2013 


根据 表 1 和 表 3 ,平缓 区 域 的 耕地 不 断 转 为 建设 用 地 ,而 
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地 类 转移 不 同 坡度 转 出 面积 The area come from different slope transfer-out 

Land type transfer 0—2° 2—6° 6—15° 15—25° 三 29? 合计 
耕地 转 林 地 Cultivated land transfer to forest land 147.38 98.87 165.05 128.37 101.36 641.03 
耕地 转 草地 Cultivated land transfer to grass land 8.07 4.22 3.52 1.05 0.16 17.01 
耕地 转 水 域 Cultivated land transfer to waters 128.56 42.41 8.61 0.92 13 180.62 
耕地 转 建设 用 地 Cultivated land transfer to build-up land 658.47 266.87 32.22 3.82 1.00 962.39 
耕地 转 未 利用 地 Cultivated land transfer to unused land 0.01 0.00 0.01 0.00 0.03 0.06 
林地 转 耕 地 Forest land transfer to cultivated land 80.60 69.32 66.37 15.26 3.51 235.07 
林地 转 草地 Forest land transfer to grass land 7.33 4.54 11.83 9.89 8.50 42.10 
林地 转 水 域 Forest land tansfer to waters 8.49 6.83 6.72 1.64 0.42 24.09 
林地 转 建 筑 用 地 Forest land transfer to build-up land 12.45 9.55 10.09 4.69 2.60 39.36 
草地 转 耕 地 Grass land transfer to cultivated land 70.84 68.56 49.49 6.65 1.30 196.85 
草地 转 林 地 Grass land transfer to forest land 29.05 24.80 48.01 23.01 9.06 133.93 
草地 转 水 域 Grass land transfer to waters 5.27 1.46 0.10 0.00 0.00 6.83 
草地 转 建 筑 用 地 Grass land transfer to build-up land 9.31 9.27 S83 2.60 1.41 28.24 
水 域 转 耕 地 Waters transfer to cultivated land 27.12 13.61 2.77 0.23 0.02 43.75 
水 域 转 林地 Waters transfer to forest land 4.34 3.37 1.94 Q.55 0.14 10.34 
水 域 转 草 地 Waters transfer to grass land 2.92 1.28 0.27 0.02 0.00 4.50 
水 域 转 建筑 用 地 Waters transfer to build-up land 2.80 1.67 0.39 0.08 0.02 4.96 
建筑 用 地 转 耕 地 Build-up land transfer to cultivated land 2.24 9.55 6.21 1.66 1.02 20.69 
建筑 用 地 转 林地 Build-up land transfer to forest land 1.15 2.89 1.37 0.31 0.12 5.83 
建筑 用 地 转 草 地 Build-up land transfer to grass land 0.27 0.32 0.14 0.02 0.00 0.75 
建筑 用 地 转 水 域 Build-up land transfer to waters 2.36 1.43 0.37 0.02 0.01 4.18 
未 利用 地 转 水 域 Unused land transfer to waters 0.00 0.00 0.00 0.00 0.00 0.00 


草地 、 林 地 又 补给 耕地 ,造成 耕地 不 断 向 山区 扩 
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张 现象 。 该 现象 在 2005 年 之 前 尤为 明显 ,2010 年 左右 出 现 耕 地 转 林 地 、 耕 地 转 水 域 现象 ,林地 水域 面积 
所 增加 ,说明 为 了 维护 生态 效应 开始 退耕 还 林 还 草 这 样 一 个 人 类 生态 修复 重建 过 程 。 
2.3 ”景观 格局 时 空 演 变 分 析 

景观 指数 高 度 浓缩 景观 格局 信息 ,能 够 反映 景观 结构 组 成 和 空间 配置 等 方面 特征 '" ,是 定量 描述 景观 格 
局 演变 及 生态 过 程 影响 的 重要 方法 '21 。 由 于 对 斑 块 的 分 析 常常 是 以 地 类 为 单位 ,分 析 某 类 型 景观 的 斑 块 特 
征 , 因 此 ,本 文选 取 类 别 层次 的 地 类 形状 指标 (FARC) 和 斑 块 聚合 度 指 标 (ALI) , 景观 层次 的 草 延 度 
(CONTAG ) Shannon 景观 多 样 性 (SHDI) .Shannon 均匀 度 (SHEI) 和 散布 与 并 列 指 标 (IJI) 进行 研究 区 23 年 景 
观 格局 变化 分 析 。 

从 图 2 各 地 类 FARC 指数 来 看 ,历年 耕地 分 维 数 最 高 ,表明 耕地 斑 块 形状 不 规则 程度 最 高 ,原因 是 :耕地 
分 布 最 广 , 且 转 出 量 最 多 ,导致 破碎 程度 高 ,但 其 面积 大 (67% 以 上 ) ,所 以 仍 显示 出 高 的 聚合 度 。 从 时 间 演 变 
来 看 ,耕地 草地、 建筑 用 地 、 林 地 斑 块 形状 随时 间 推 移 愈加 不 规则 ,水 域 和 未 利用 地 形状 则 呈 相 反 发 展 趋势 ; 
林地 在 2010 年 后 聚合 度 下 降 .而 水 域 和 草地 呈 高 的 聚集 状态 。 

从 表 4 景观 尺度 指标 来 看 ,研究 区 CONTAG 呈 逐 年 下 降 趋 势 , 由 1990 年 的 68.10 降 为 2013 年 的 67.59。 
SHDI SHEI 和 IJ 均 逐 年 上 升 , 且 JI 上 升幅 度 高 于 SHDI,SHDI 高 于 SHEI, 这 是 水 域 .草地 等 地 类 形状 指标 趋 
于 规则 .聚集 程度 上 升 共 同 作用 的 结果 。SHDI 和 SHEI 值 上 升 ,表明 23 年 间 地 类 破碎 现象 一 直 在 延续 ,优势 
地 类 不 断 减弱 ,地 类 分 布 复 杂 化 。CONTAG 逐年 下 降 , 即 为 优势 度 下 降 所 导致 的 景观 间 连 接 性 和 通达 性 的 
下 降 。 
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号 

Eee 1995FRAC AM 
三 了 14 me 2000FRAC AM 
名 ] 95 mm 2005FRAC AM 
又 号 首 mm 2010FRAC AM 
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区 六 长 mm 2013FRAC AM 
秘 当 」] 90 法 
兴 旦 赋 咏 ”一 @ 一 1990AI 
家 中 12 《<2 
RR 妃 加 ”一 @ 一 1995AI 
E | 85 4 ~。 2000AI 
i | 
便 宇 1 | 一 e@_ 2005AI 

oo 

隔 | | | —®— 2010AI 

1.0 | - 本 80 —® 2013AI 
耕地 ”建筑 用 地 林地 水 域 草地 荒漠 
土地 利用 类 型 Land use type 


图 2 淮河 上 游 1990 一 2013 年 不 同 土地 利用 类 型 FRAC_AM 和 AI 指数 变化 分 析 
Fig.2 Changing analysis different land use FRAC AM and AI index at the upper reaches of the Huaihe River from 1990 to 2013 


表 4 淮河 上 游 1990 一 2013 年 景观 格局 指数 表 
Table 4 The landscape pattern index table of the upper reaches of the Huaihe River from 1990 to 2013 


. re jn 1990 1995 2000 2005 2010 2013 
CONTAG 68.1046 68.0436 68.022 67.7877 67.8589 67.5849 
IJI 49.8598 49.8612 49.891 50.1076 51.2142 51.2435 
SHDI 1.5616 1.563 1.5665 1.5789 1.6173 1.6329 
SHEI 0.5512 0.5514 0.5529 0.5573 0.5595 0.565 


CONTAG :蔓延 度 ， contagion ; JI: 散布 与 并 列 指数 , interspersion juxtaposition index; SHDI: 香农 多 样 性 , shannon’s diversity index; SHEI: 香农 均匀 


度 ,shannon’s evenness index 
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3 ESV 变化 分 析 


生态 系统 的 异 质 性 .生态 系统 服务 的 复杂 性 ,造成 ESV 测量 困难 性 55] 。Gret-Regameyl23] 等 认为 ESV 研 
究 不 是 给 区 域 生 态 环境 加 上 一 个 准确 的 、 绝 对 的 价格 标签 , 而 是 探讨 单位 自然 资产 增加 或 减少 时 其 服务 价值 
的 边际 变化 ,进而 探讨 一 种 合理 的 土地 利用 格局 。 已 有 研究 表明 ,边际 变化 分 析 更 适合 于 研究 ESV 
变化 '**”3]。 
3.1 ESV 时 间 变 化 特征 分 析 

参照 河南 省 1990 一 2013 年 统计 年 鉴 ,进行 主要 粮食 作物 播种 面积 及 产值 的 计算 。 在 此 基础 上 ,根据 各 年 
份 CPI 和 式 (3) 一 (6) ,得 到 剔除 通货 膨胀 部 分 的 单位 面积 价值 ,结合 各 生态 系统 面积 数据 和 式 (1) ,计算 研究 
区 各 年 份 ESV( 表 5) 。 


表 5 淮河 上 游 1990 一 2013 年 ESV 时 间 变 化 表 
Table S The time variation of ESV of the upper reaches of the Huaihe River from 1990 to 2013 


生态 系统 服务 各 年 份 ESV Each year ESVZ108 元 

Ecosystem services 1990 1995 2000 2005 2010 2013 

供给 服务 食物 生产 27.36 22.39 22.81 19.27 31.31 38.80 
Provisioning service 原料 生产 12.94 9.59 10.66 8.98 14.78 18.34 
水 资源 供给 -6.17 -10.20 一 4.94 -3.40 -5.12 一 6.42 

调节 服务 气体 调节 34.93 27.89 28.83 24.31 40.77 50.67 
Rrgulating service 气候 调节 54.92 42.67 44.99 38.01 66.38 82.67 
净化 环境 20.23 15.20 16.89 14.70 25.51 31.75 

水 文 调节 126.26 95.25 111.08 103.24 182.00 225.49 

支持 服务 土壤 保持 38.83 28.39 31.86 26.85 44.69 55.50 
Support service 维持 养分 循环 4.90 3.95 4.05 3.42 5.65 7.01 
生物 多 样 性 23.73 17.67 19.51 16.52 27.81 34.73 

文化 服务 Cultural service ”美学 景观 11.42 8.37 9.43 8.03 13.40 16.74 
合计 Total 349.34 261.16 295.18 259.92 447.18 555.27 


由 表 5 可 知 ,1995 一 2005 年 ESV 波动 下 降 ,2005 年 后 呈 明 显 上 升 趋势 。11 种 生态 系统 服务 功能 亦 呈 现 
上 述 年 际 变化 特征 。 其 中 调节 服务 功能 增值 最 多 ,其 次 为 支持 服务 功能 ,分 别 增加 154.24 亿 元 和 29.79 亿 元 。 
结合 表 6,6 类 ESV 在 1995 一 2005 年 间 均 有 不 同 程度 下 降 ,2005 年 之 后 明显 上 升 。 林 地 水域 和 农田 是 ESV 
的 主要 构成 部 分 ,23 年 间 ,水 域 增值 最 大 ,增值 97.55 亿 元 ,其 次 为 林地 农田 ,分 别 增加 72.70 亿 元 和 39.12 亿 
元 ,湿地 为 唯一 负增长 项 ,减少 4.64 亿 元 。 


表 6 淮河 上 游 1990 一 2013 年 各 生态 系统 ESV 表 
Table 6 The ESV of the upper reaches of the Huaihe River from 1990 to 2013 


生态 系统 各 年 份 ESV Each year ESV/103 元 

Ecosystem 1990 1995 2000 2005 2010 2013 

农田 Farmland 101.50 80.36 84.45 71.09 113.62 140.62 
林地 Forest land 127.25 105.55 105.54 88.99 160.39 199.95 
草地 Grass land 14.70 4.23 9.42 7.67 12.77 15.90 
水 域 Waters 89.42 62.90 82.45 81.94 151.66 186.97 
湿地 Wet land 16.47 8.12 13.30 10.22 8.75 11.83 
荒地 Waste land 0.00 0.00 0.00 0.00 0.00 0.00 

合计 Total 349.34 261.16 295.18 259.92 447.18 555.27 


根据 表 1 和 表 6, 水 域 和 林地 所 占 区 域 面 积 比例 分 别 为 3.26% .17.62% ,ESV 贡献 量 却 占 有 33.67% 和 
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36.00% ,表明 水 域 和 林地 是 引起 ESV 提升 的 主要 土地 利用 类 型 。 耕 地 、 草 地 面积 的 减少 ,其 ESV 却 在 增加 ， 
这 一 现象 的 产生 在 于 研究 区 主要 粮食 单位 产值 的 增加 ,致使 各 生态 系统 单位 面积 价值 增加 ,抵消 了 面积 减少 
带 来 的 负面 效应 。 
3.2 生态 系统 服务 损益 分 析 

基于 式 2 和 1990 一 2013 年 土地 转移 矩阵 进行 ESV 损益 计算 ,得 到 表 7。 由 表 7 可 知 23 年 间 ,ESYV 增加 
205.93 亿 元 ,主要 由 农田 转 为 林地 农田 转 为 水 域 带 来 的 生态 系统 服务 增值 构成 ;所 有 地 类 转 水 域 均 为 正 流 
向 ,增值 178.52 亿 元 。ESV 损失 主要 由 林地 草地 水域 转 农田 造成 ,林地 转 农 田 损 失 最 多 , 达 138.05 亿 元 。 
所 有 地 类 转 建筑 用 地 均 呈 负 流 向 ,损失 80.34 亿 元 ,该 部 分 损失 主要 由 农田 转 建设 用 地 构成 。 整 体 来 看 ， 
1990 一 2013 年 ESV 呈 收 益 大 于 损失 情况 ,生态 增值 主要 来 源 于 林地 和 水 域 生 态 价值 贡 献 。 


表 7 淮河 上 游 1990 一 2013 年 ESV 损益 矩阵 


Table 7 The ecosystem services profit or loss maxtrix values of the upper reaches of the Huaihe River from 1990 to 2013 


土地 利用 类 型 农田 林地 草地 水 域 建筑 用 地 总 计 
Land use type Farmland Forest land Grass land Waters Build-up land Total 
农田 Farmland/108 元 0.00 127.97 26.16 82.22 -38.32 198.03 
林地 Forest land/103 元 -73.91 0.00 -16.51 36.64 -35.35 —89.14 
草地 Grass land/108 元 -78.85 48.35 0.00 24.77 —0.55 -6.27 
水 域 Waters/10* 元 -53.22 —17.15 —5.24 0.00 -6.13 —81.74 
建筑 用 地 Build-up land/10* 元 67.92 51.99 30.25 34.89 0.00 185.06 
总 计 Total/105 元 —138.05 211.16 34.66 178.52 —80.34 205.93 


3.3 ”ESV 空间 变化 特征 分 析 

从 研究 区 地 形 分 布 来 看 , 沿 西北 -东南 向 ,由 南 至 北 , 依 次 呈 山 地 .丘陵 .平原 三 个 梯 段 ,图 3 各 县 市 ESV 亦 
呈现 梯 段 变化 规律 :山区 面积 居多 的 泌 阳 县 、 桐 柏 县 、 信 阳 市 . 罗 山 县 .新 县 和 商城 县 6 个 时 期 的 ESV 普遍 最 
高 ,丘陵 面积 居多 的 确 山 县 .光山 县 江川 县 和 固 始 县 6 个 时 期 的 ESV 较 高 ,平原 区 域 的 县 市 ESV 较 低 。 其 
中 , 汝 南 县 因 其 水 资源 丰富 ,ESV 偏 高 于 其 他 平原 区 域 县 市 。 各 县 市 1995 年 和 2005 年 ESV 均 呈 明显 低谷 
区 ,2005 年 之 后 增值 迅速 ,2010 年 之 后 ,山区 县 市 仍 呈 高 速 增长 趋势 ,其 他 区 域 县 市 2013 年 ESV 基本 与 2010 
年 相 平 。 驻 马 店 市 是 唯一 的 负增长 区 域 , 究 其 原因 ,城市 建设 迅速 ,地 类 多 样 性 下 降 严 重 。 

结合 表 7 和 图 3 ,23 年 间 生态 系统 服务 总 值 增 加 ,增值 主要 来 源 于 林地 与 水 域 生态 价值 贡献 ,而 林地 与 水 
域 的 分 布 多 依赖 于 地 形 因素 , 因 此 出 现 生态 系统 服务 总 值 增加 ,但 又 空间 分 布 不 平衡 现象 , 且 耕 地 不 断 被 转移 
为 建筑 用 地 所 造成 的 负面 效应 已 经 呈现 ,因此 未 来 应 更 加 重视 自然 资源 保护 和 生态 服务 价值 的 空间 均衡 化 ， 
使 得 人 人 得 以 享受 生态 环境 改善 所 带 来 的 增值 服务 。 


4 ESV 与 土地 类 型 变化 .土地 景观 变化 的 响应 关系 分 析 


4.1 土地 利用 变化 与 ESV 变化 的 相关 性 分 析 

采用 研究 区 1990 年 .1995 年 .2000 年 .2005 年 .2010 年 和 2013 年 6 个 年 份 的 一 级 地 类 面积 与 相应 ESV 
进行 相关 性 分 析 ,相关 系数 为 0.35 ,在 0.10 水 平 上 达到 显著 。 表 明 经 过 CPI 单位 价值 系数 修正 后 的 ESV 与 6 
大 地 类 面积 存在 一 定 的 相关 性 。 那 么 6 大 地 类 的 空间 分 布 情况 与 其 ESV 存在 怎样 的 相关 性 ?以 各 县 市 为 统 
计 单 元 (图 3) ,以 1990 年 和 2013 年 的 土地 利用 面积 变化 量 与 ESV 变化 量 进行 相关 性 分 析 , 其 中 ,变化 量 为 
前 .后 两 个 时 期 的 差 值 ,两 组 变化 量 的 相关 系数 为 0.38 ,在 0.01 水 平 上 达到 显著 ,相关 程度 明显 高 于 两 个 总 量 
的 相关 性 ,表明 ESV 在 随 面积 变化 影响 的 基础 上 ,主要 受 地 类 面积 的 空间 分 布 影响 。 那 么 地 类 的 时 空 演 变 过 
程 是 否 与 ESV 存在 一 定 联系 ? 仍 以 县 市 为 统计 单元 ,以 1990 年 和 2013 年 的 土地 利用 变化 率 与 ESV 变化 率 
进行 相关 性 分 析 , 变 化 率 参照 文献 "1( 式 (3-3) ) 获得 ,相关 系数 为 0.28 ,在 0.05 水 平 上 达到 显著 。 该 值 略 低 
于 ESV 与 空间 分 布 变化 量 相 关 性 ,表明 地 类 变化 的 时 间 效 应 对 ESV 的 影响 程度 要 低 于 地 类 的 空间 位 置 变化 。 
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3 淮河 上 游 1990 一 2013 年 各 县 市 ESV 空间 分 布 图 
Fig.3 The spatial distribution of ESV in each counties of the upper reaches of the Huaihe River from 1990 to 2013 


但 研究 区 地 形 差 异 大 ,处 于 南北 气候 交接 带 ,自古 以 来 是 我 国 重 要 的 粮食 棉花、 油料 生产 基地 ,农业 生 
产 水 平和 土地 利用 率 较 高 ,各 地 类 转移 频繁 ,不 同 地 类 的 景观 指数 与 ESV 又 表现 为 怎样 的 相关 性 ? 基于 此 ， 
以 下 从 土地 景观 格局 角度 ,分 析 不 同 地 类 景观 指标 与 ESV 的 响应 关系 。 
4.2 景观 格局 变化 与 ESV 变化 的 相关 性 分 析 

利用 Fragstats 4.2 软件 进行 斑 块 个 数 ( NP ) 、 斑 块 平均 面积 ( AREA_MN) .最 大 斑 块 所 占 景观 面积 的 比例 
(LPI) 形状 (SHAPE_MN) .形状 维 数 ( FRAC_MN) .聚集 户 (AI)6 个 景观 指数 的 计算 ,分 别 得 到 出 6 期 农田 、 
林地 草地 水域 .湿地 和 荒地 6 个 地 类 的 景观 指数 。 在 此 基础 上 ,对 6 个 地 类 的 每 个 指数 与 各 期 ESV 进行 相 
关 性 分 析 。 得 出 下 表 : 


表 8 淮河 上 游 1990 一 2013 年 景观 格局 指数 与 ESV 相关 性 
Table 8 The relevance of landscape pattern index and ESV at the upper reaches of the Huaihe River from 1990 to 2013 


景观 指数 斑 块 数量 面积 平均 。 最 大 斑 块 指数 ” 斑 块 形状 指数 分 形 维 聚集 度 
Landscape index NP AREA_MN LPI SHAPE_MN FRAC_MN AI 
相关 系数 Correlation coefficient .73"* 0.35” 0.33* 一 0.09 一 0.34 0.56** 


**0.01 水 平 上 显著 相关 ; * 0.05 水 平 上 达到 显著 


表 8 表明 ESV 与 景观 丰 度 和 地 类 聚集 度 存在 明显 相关 性 ,与 地 类 闹 块 大 小 相关 性 偏 弱 ,与 各 斑 块 形状 复 
杂 程 度 呈 较 强 负 相 关 性 ,说 明了 景观 破碎 化 对 ESV 的 影响 程度 。 
4.3 ESYV 时 空 变化 的 影响 因素 分 析 
en 研究 表明 ,土地 利用 变化 影响 生态 系统 服务 ,但 对 影响 的 过 程 和 机 理 认识 仍然 不 足 '”1。 相 关 分 析 是 
研究 空间 关系 问题 ,分 析 空 间 格局 的 重要 方法 和 有 效 手段 。 运 用 相关 系数 量化 土地 利用 指标 与 ESV 的 
2 ESV 动态 变化 规律 ,进而 探讨 构建 更 为 合理 的 生态 系统 空间 格局 。 研 究 表明 :(1)ESYV 与 地 
类 空间 分 布 相关 性 显著 , 表 3 和 图 3 的 空间 变化 一 致 性 亦 验 证 了 该 结论 ,也 与 胡 和 兵 等 .| 和 欧阳 志 云 等 引 研 
究 结 果 相 一 致 。 均 表现 为 受 地 形 地 貌 影响 的 水 热 条 件 好 、 利 于 植被 生长 的 区 域 ,其 生态 系统 的 单位 面积 价值 
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量 较 高 ,而 往往 这 些 区 域 景观 多 样 性 和 生态 完整 性 也 较 好 ,因此 ESV 与 NP 和 AI 亦 呈 明显 相关 性 。(2) ESV 
与 地 类 年 际 变 化 相关 性 偏 高 。 从 表 1 可 知 ,1995 一 2005 年 地 类 转移 最 为 频繁 ,大 面积 耕地 转 为 建筑 用 地 ,而 
林地 和 草地 转 出 补给 耕地 ,2010 年 后 ,部 分 耕地 转 出 为 林地 和 草地 。 而 ESV 在 相应 年 份 分 别 出 现 波 动 下 降 ， 
和 明显 上 升 趋势 ( 表 6) ,两 者 时 间 节 点 变化 一 致 。 虽 然 林地 与 水 域 面积 比例 小 ( 表 1) , 受 ESV 当量 影响 ,两 者 
对 区 域 ESV 总 值 贡 献 比 重大 ,成 为 ESV 的 关键 “热点 ”区 '*”i。 也 验证 了 ,对 于 不 同类 型 的 生态 系统 ,面积 变 
化 与 其 ESV 相关 性 偏 弱 这 一 结论 。(3)ESV 与 AREA_MN 和 LPI 相关 性 偏 高 ,与 FRAC_MN 存在 负 相 关 ,与 
孝 芳 华 等 '50 研究 结果 相 一致 。 这 是 因为 ,23 年 间 景观 破碎 化 在 逐年 加 重 ,景观 优势 度 逐 年 降低 ,这 一 现象 在 
1995 年 左右 表现 最 为 明显 (图 3) ,生态 系统 “热点 ”区 域 不 断 被 破坏 ,也 说 明了 人 为 活动 是 影响 土地 景观 格局 
的 主要 因素 ,而 城镇 扩张 引发 生态 用 地 的 减少 往往 是 导致 ESV 总 值 损失 的 主要 原因 5 。(4) 研究 地 类 结构 
和 格局 时 空 演变 特征 , 细 化 地 类 结构 .格局 指标 ,并 按 各 指标 与 ESV 的 相关 程度 进行 梳理 ,分 析 自 然 和 人 为 因 
素 在 ESV 过 程 变化 中 的 作用 机 理 , 较 好 的 克服 了 已 往 评 价 尺度 过 大 ,评价 内 容 过 泛 '” ,结论 针对 性 不 强 的 问 
题 。 成 为 探究 土地 利用 变化 下 生态 系统 服务 响应 特征 研究 的 一 种 有 益 尝试 。 今 后 仍 需 加 强 土地 利用 、 景 观 空 
间 格 局 .CPI 对 最 终 的 ESV 影响 定量 判定 研究 ,以 深入 揭示 土地 利用 变化 对 生态 系统 服务 的 影响 过 程 。 


5 结论 


ESV 与 地 类 的 的 空间 分 布 属性 相关 性 最 为 显著 ,与 地 类 的 变化 速度 相关 偏 弱 ;ESV 与 土地 景观 异 质 性 和 
破碎 性 相关 性 显著 ,对 景观 的 形状 特征 相关 性 偏 弱 。 

采用 生态 系统 价值 损益 分 析 方 法 ,使 价值 变化 更 多 关注 在 自然 资源 的 过 渡 性 损失 方面 ,弥补 了 现 有 生态 
系统 服务 在 过 程 变 化 评价 方面 的 欠缺 ,更 利于 评价 结果 在 土地 利用 实际 管理 中 的 运用 。 

采用 CPI 指数 修订 ESV 系数 ,能 够 消除 时 间 序 列 内 各 年 份 物价 上 涨 带 来 的 价值 变化 ,避免 了 价格 影响 下 
的 ESV 研究 结果 不 尽 一 致 问 题 ,使 得 计算 结果 更 贴近 实际 生态 系统 价值 ,提高 了 多 时 间 序 列 生态 系统 服务 评 
价 的 准确 性 和 可 比 性 。 
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